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Abstract 

We have selected 42 elliptical galaxies from the literature and estimated their ve- 
locity dispersions at the effective radius (ctrJ and at 0.54 effective radii (o"o.54R e ). 
We find by a dynamical analysis that the normalized velocity dispersion of the dark 
halo of an elliptical galaxy <tdm is roughly an a multiplied by a constant, which is 
almost independent of the core radius or the anisotropy parameter of each galaxy. 
Our sample analysis suggests that <r^ )M lies in the range 178-198 km s -1 . The power 
law relation we find between the luminosity and the dark matter velocity dispersion 
measured in this way is (L/L*) = (odm/odm) 7 ' where 7 is between 2-3. These 
results are of interest for strong gravitational lensing statistics studies. 

In order to determine the value of crJj) M , we calculate Mg T ^ in the same Bt (0) 
band in which o"J) M has been estimated. We select 131 elliptical galaxies as a com- 
plete sample set with apparent magnitudes Bt (0) between 9.26 and 12.19. We find 
that the luminosity function is well fitted to the Schechter form, with parameters 
Mg T j j = -19.66 + 5- log 10 /i± 0.30, a = 0.15 ±0.55, and the normalization constant 
4>* = (1.34 ± 0.30) x 10~ 3 /i 3 Mpc~ 3 , with the Hubble constant H Q = 100 h km 
s _1 Mpc -1 . This normalization implies that morphology type E galaxies make up 
(10.8 ± 1.2) per cent of all galaxies. 
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1 INTRODUCTION 



The search for the dark matter plays a central role in modern astrophysics and cosmology. A 
key factor in constraining the cosmic abundance of dark matter is to measure rotation curves 
or velocity dispersions of galaxies. Because spiral galaxy rotation curves provide sensitive 
constraints on dark matter, it is common to measure velocity dispersions in these systems. 
However, there are fewer available data sets of velocity dispersions for elliptical galaxies. 

At the same time, the characteristic velocity dispersion of elliptical galaxies is a very impor- 
tant parameter in the statistical study of strong gravitational lensing, because the optical 
depth for lensing by galaxies is proportional to the fourth power of the velocity dispersion 
(Turner, Ostriker & Gott 1984). The value of the normalized dark matter velocity dispersion 
for luminous elliptical galaxies that should be used in the strong gravitational lensing anal- 
ysis is still under some debate. For example, in 1984 Turner et al. (1984) suggested that one 
should use directly observed velocity dispersions, and multiply them by ^3/2 to determine 
dark matter velocity dispersions, as suggested by the virial theorem. However, Kochanek 
(1993,1994) has argued using dynamical models that one should simply use the measured 
velocity dispersions as the dark matter velocity dispersions for the lensing calculation, with- 
out any ^3/2 correction. However, even in this case, one must still estimate the characteristic 
dark matter velocity dispersion for ellipticals, and this depends on one's assumptions about 
the elliptical galaxy luminosity function (Cheng & Krauss 2000). Moreover, both Kochanek 
(1993,1994) and Turner et al. (1984) focussed on the singular isothermal sphere (SIS) model 
in their gravitational lensing analysis. This model is clearly a very rough approximation to 
actual galaxies. 

In this paper, our goal is to re-examine both the appropriate luminosity function and veloc- 
ity dispersions for use in lensing statistics analyses. In our analysis we consider a finite core 
model, and consider the self-consistent determination of the velocity dispersion and lumi- 
nosity function for ellipticals. This analysis is relevant even if the core radius is quite small, 
or zero. To be consistent we calculate the luminosity function in the same band in which we 
estimate velocity dispersions. 

We choose the Hubble constant H Q = 100 h km s _1 Mpc -1 throughout. We investigate the 
luminosity function in Section 2. In section 3 we review our investigation of the theoretical 
determination of the dark matter velocity dispersion in galaxies with core radii using dy- 
namical models with luminosity profiles discussed by Jaffe (1983) and Hernquist (1990). In 
Sections 3.1 and 3.2, we will describe our main results appropriate for use in lensing studies, 
including normalized velocity dispersions and their power law relations with the luminosity. 
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2 THE LUMINOSITY FUNCTION 



The appropriate value of M* in a Schechter formulation (1976) of the luminosity function 
will depend on which magnitude system one is measuring, so that there is no unique value 
of M*. In this paper, we will investigate the shape of the Schechter luminosity function for 
a sample of elliptical galaxies in the Bt (0) band, as defined in de Vaucouleurs et al. (1991) 
(hereafter RC3). In order to calculate the luminosity of each galaxy in the Bt (0) band, we 
choose the distance of each galaxy from Faber et al. (1989) (hereafter 7s), and the apparent 
magnitude B T (0) of each galaxy from RC3. We find 415 galaxies in this set. However, 
we need to obtain a complete magnitude limited sample set to calculate the luminosity 
function (Loveday et al. 1992). Based on the discussion in 7s, the elliptical galaxies in 7s are 
complete to apparent magnitude 11.89. However, the RC3 catalogue was published after the 
appearance of 7s. Therefore, we consider galaxies with the apparent magnitudes given in 7s 
less than or equal to 11.89, and then check out their corresponding apparent magnitudes 
in RC3. We find the dimmest (or the largest number) apparent magnitude among these 
galaxies is 12.19. Thus, we use 12.19 as the cutoff apparent magnitude in RC3, and select 
those galaxies in 7s whose corresponding apparent magnitudes in RC3 are less than 12.19. 
There are 131 galaxies in the final sample set. We then determined the luminosity function 
of this sample as described below. 



2.1 STY parametric maximum likelihood method 



We use here STY maximum likelihood method (Sandage, Tammann & Yahil 1979). There 
are some advantages to use this method, as discussed in Efstathiou, Ellis & Peterson (1988) 
(hereafter EEP). Consider the following probability 

<t>{ L i) m 
Pi oc — f m (1) 



with Schechter (1976) function 



^(^) = ^(^) a exp(-^) (2) 



This is the probability to find a galaxy at distance di with luminosity Lj in the magnitude- 
limited sample set. L max (rfj) and L min (rfj) are calculated by using the distance di (in unit 
of h~ x Mpc) and the brightest and the dimmest apparent magnitudes among the whole 
magnitude-limited sample set. For example, the brightest apparent magnitude in our sample 
set is 9.26, so L max (d^ /L* will be equal to dj ■ io 10 +°- 4 ( M *- m — ), with m max = 9.26. The 

Af 

likelihood function is C = \\pi, and Af is the total number of galaxies in the sample set. 

i=l 
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Fig. 1. Parameters, a and M Bt ^ — 5 • log 10 h, are fitted to the Schechter function by using STY 
maximum likeiihood method. The best fit has been marked '+.' The solid curve is the 1 a error 
contour, and the dash curve is the 95 per cent confidence level contour. The statistical method is 
discussed in the text. 

We maximize the likelihood function C with two parameters a and L* (or M*). We do not 
additionally correct for a Malmquist bias because the average dispersion between different 
B T magnitudes in RC3 is small, and also because Faber et al. (1989) already corrected for 
the Malmquist bias when they estimated the final distances of the galaxies in 7s. From the 
maximum likelihood analysis, we get the parameters = -19.66 + 5-log 10 /i lo33 an d 

a = 0.15 toil (1 o error) (Figure 1). 

Because of the use of two different samples to attempt to obtain a complete magnitude 
limited survey, one might be concerned that different choices of magnitude cutoff would 
produce different maximum likelihood estimates for the luminosity function parameters. We 
explored a variety of combinations in this regard, and found consistent results within the 
uncertainties quoted above. For example, if we consider galaxies whose apparent magnitudes 
are below the apparent magnitude threshold in 7s, then find the corresponding apparent 
magnitudes in RC3, and neglect one apparently anomalous faint galaxy in 7s with apparent 
magnitude 11.00 in RC3, then the next galaxy which is below the threshold in 7s has apparent 
magnitude 11.83 in RC3. If we use 11.82 as the cutoff apparent magnitude, then we will have 
82 galaxies in the sample set, and then we get M^™ = -19.78 and a = 0.22. Choosing a 
fainter cutoff, 11.98, we find Mg T ^ = -19.74 and a = 0.00 with 107 galaxies in the sample 
set. As claimed, these are within the uncertainties in our quoted value above. 
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2.2 Step-wise maximum likelihood method 



In order to examine the goodness of fit of the Schechter function to the data we utilize a 
Step-wise maximum likelihood method to fit the data (EEP). This method is basically to 
replace the integral in equation (1) with finite sum. Following presentation in EEP, we define 



where 



N 



Y^De (Li-L k ) 



4>'k 



1=1 



N 



Th [L k , L m i n (di) , L m3iX (di)] 



5 Np 

<P'j ■ Th [Lj, L min (di) , L max (di)] 

3=1 



(3) 



<t>' k (L) 
De(Li - L k ) 
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Here N is the total number of galaxies in the magnitude-limited sample set, and N p is the 

total number of bins with bin width AL. The index k runs from 1, 2, 3, to iV p . Here the 

step function Th[L k , L min (di) , L max (di)], is different from that in EEP, because the upper 
limit of the integral in equation (1) is L max (di) rather than infinity. We simply use the results 
obtained in the previous subsection, put them in equation (3) to compare the values (along 
with errors) for the discrete points derived in the second method to the Schechter function 
form (Figure 2). We have to take some care to normalize the values of (f)' k , if we want to 
compare these two maximum likelihood fits. The uncertainty of each dot is obtained using 
equation (3), using straightforward error propagation (Bevington & Robinson 1992). 

We can see from the figure that these two methods are consistent with each other for the 
Schechter parameters from the previous subsection, implying that the Schechter function is 
a good fit to the data. This remained true for other choices of bin widths, AL/L* = 1 and 
AL/L* = 0.1 as well. 
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Fig. 2. Comparison between the step-wise maximum likelihood method (discrete dots) and the STY 
method (solid curve). We use results from Section 2.1, and bin width AL/L* = 0.5 for the step-wise 
method in this plot. 

2.3 The normalization constant 



We calculate the normalization constant 0* using standard techniques (EEP or Loveday et 
al. (1992)). This involves calculating the average space density of galaxies in the luminosity 
range being probed, n, as follows: 



n 



N 

£ 

i=i 



w (di 



jj max dx4irx 2 S (x) w (x) 



(4) 



where 



rinin( Lmaic (d» ) ,L 2 ) , ^ rfL 

S%4>{L)dL 
1 

1 + 4nnJ 3 (r c ) 5 (d*) ' 



S(di) 
w (di) 



Js (r c ) 



I c 

/ r 2 e (r) 



iV is the total number of galaxies in the sample set, di is the distance of each galaxy, and d m i n 
and <i max are the nearest and the furthest distances among galaxies in the sample set. We have 
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c^min = 6.95 h^ 1 Mpc, and d max = 57.79 h~ l Mpc. L\ and L2 are related to the selection func- 
tion S(di) and the weighting function w(di), but these can, with sufficient accuracy for our 
purposes be chosen to be L min (dj) and L max (dj), the dimmest and the brightest luminosities 
from our sample set. Finally, £(r) is the two-point galaxy correlation function with r c 20/i _1 
Mpc, and we use 47rJ 3 ps 10,000/i" 3 Mpc 3 here (Efstathiou et al. 1988,Loveday et al. 1992). 

After applying the Schechter function in the selection function S(di), we find n = 9.65 x 
10 _4 /i 3 Mpc -3 from equation (4). The normalization constant is equal to 



III <f> ( L ) dL 



(5) 



So we get 4>* = 1.34 x 10 3 h 3 Mpc 3 . The uncertainty of n is approximated by (Loveday et al. 1992) 



Sn 



n 



\ l? max wSdv 



(6) 



We thus find 5n/n 0.19. This value is larger than that in Loveday et al. (1992), because 
#n/n is inversely proportional to the square root of number of galaxies in the sample set. 
There are 131 galaxies in our sample set, but there were about 1600 galaxies in Loveday et 
al. (1992). We thus find the uncertainty of 0* by standard methods of error propagation, 
and hence obtain our final result, (jf = (1.34 ± 0.30) x 10~ 3 h 3 Mpc" 3 . 



3 MODELS OF VELOCITY DISPERSION 



We have already mentioned the importance of using a self-consistent analysis of dark matter 
velocity dispersions when attempting to model galaxies as lenses (Cheng & Krauss 2000). 
Here we describe such an analysis for various galaxy models. The key point is that the 
determination of o"dm is a sensitive function of the galaxy anisotropy parameters and core 
radius, and one must consider this sensitivity when determining o"dm o n the basis of various 
observations. A theoretical dynamical analysis suggests that measuring line of sight velocity 
dispersions at the de Vaucouleurs effective radius R e (1948), can provide a relatively robust 
estimate of the underlying <7dm, and also allow one to extract out contributions due to central 
mass concentrations in the galaxy. 

Throughout our analysis we model the mass density distribution of elliptical galaxies with 
the following form (Hinshaw & Krauss 1987): 

2 

_ °"dm 
9 ~ 2nG (r 2 + r 2 c ) 



(7) 



where <tdm (independent of radius r) is the velocity dispersion of this system (presumably 



7 



Rj R c Rj R e 




0.01 0.05 0.1 0.5 1 5 10. 0.01 0.05 0.1 0.5 1 5 10. 



Rj R c Rj R e 

Fig. 3. Jaffe (1983) luminosity profile with core radii, 0, R c /20, R e /40, and R e /60. The five curves 
in each plot, dash-dot curve, dot-dot curve, solid curve, dash-dot-dot curve, and dash-dash curve, 
represent the five different anisotropy parameters (3 = 0.4, 0.2, 0.0, -0.2, and -0.4, respectively. 

the dominant dark matter), and r c is the core radius. When r c is zero, this model reverts to a 
singular isothermal sphere (SIS) model. We should bear in mind that the velocity dispersion 
in this equation, which is assumed to be independent of radius, cannot be measured directly. 
What we can measure is the line-of-sight velocity dispersion <t 1os (R), with the projected dis- 
tance R measured from the center of the observed galaxy (i.e., R is perpendicular to the line- 
of-sight to us). This means that we need to integrate over the line-of-sight distance in order 
to work backward to infer the value of <tdm for theoretical purposes. For a singular isothermal 
sphere, <tdm is not a function of R, but <xi os (R) is, and this relationship is different when a 
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Fig. 4. Hernquist (1990) luminosity profile with core radii, 0, R e /20, R e /4Q, and R e /60. The de- 
scription of the five curves in each plot is the same as in Figure 3. 

galaxy has a finite core (although this is a secondary issue). Thus, one has to be careful how 
to use the measured velocity dispersion to derive the relevant quantity to utilize for lensing, 
namely whether it well approximates odm- To find the relationship between o"dm and o"i os (R), 
we start from Jeans equations as follows (Ogorodnikov 1965,Binney & Tremaine 1987): 

_(^) +23 _, = -,_ ( 8 ) 



9 



where 



<9$ 

dr 



4ttG 
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1 



dr' 



r' 2 p (r') 



2d 



DM 



arctan 



(9) 



by using equation (7). 

In equation (8), u(r) is the luminosity profile, a T (r) is the radial part of the velocity dispersion 
of the observed galaxy at distance r from the center of that galaxy. The other two components 
of the velocity dispersion a in spherical coordinates are <jg and a^. We define the anisotropy 
P(r) as ^ = = (1 - fit*. 

The line-of-sight velocity dispersion a\ os (R) is related to va 2 by 

/ (R) al (R) = 2f(l- (5-^ ^-^L^dr (10) 

where I(R) is the surface brightness profile, which can be calculated as follows: 

oo 
R 



By solving the differential equation (8), we cam re-write equation (10) as 

°r ( R 2 \ 1 °r fl(f> 

I { R ) „ L(R) = 2j d „(l-p^ w - w J rfx.^M- (12) 



VR 2 +u 2 



If we assume f3 is simply a constant independent of r in equations (9), (11), and (12), then 
once we know the form for the luminosity profile z/(r), we know the ratio of <yf os /a'^ )M at 
the different projected distance R. We have considered two different but similar luminosity 
profiles v(r). We adopt u(r) oc r~ 2 (r + r )~ 2 , with r = 1.31i? e , from Jaffe (1983), and 
v(r) oc r _1 (r + r )~ 3 , with r = 0.45i? c , from Hernquist (1990). We performed our numerical 
calculations by adjusting f3 = 0.4, 0.2, 0, -0.2, -0.4 and r c = 0, ^R e , j^Re, for both 
Jaffe (1983) and Hernquist (1990) luminosity profiles (Figure 3 and Figure 4). The ratio 
of the core radius to the effective radius may be estimated from Lauer (1985), and we 
find that this ratio is between 1/20 and 1/60 for class I galaxies in Lauer (1985). From 
Figure 3 and Figure 4, we clearly see that of os /o" DM is strongly sensitive to the anisotropy 
when the projected radius is less than 0.1i? e (whether or not the core radius is non-zero, 
a factor to which the result is also sensitive). This is consistent with observational data 
(Bender, Saglia & Gerhard 1994). It can thus be dangerous to simply consider the line-of- 
sight velocity dispersion within 0.1i? e as <7dm- Even if the core radius is effectively very small, 
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as numerous observations of ellipticals suggest, the sensitivity to the anisotropy parameter 
implies one should be careful when interpreting meassured values of the velocity dispersion. 

However, if we are careful to look at all the 40 curves in Figure 3 and Figure 4, we can find 
the following inequalities: 

1.16 < < 1.27 ati? = 0.4i? e 

1.20 < < 1.30 at R = 0.54i? e 

1.24 < -2%r < 1.37 a tR = R c 

These indicate that the intrinsic scatter of o"dm will be less than 10 per cent if we can 
measure the line-of-sight velocity dispersion at R c or half R e . Then, we can multiply this 
velocity dispersion by the average value (e.g., 1.31 for R = R c , or 1.25 for R = 0.54i? e ), in 
order to get <tdm- We can see that this argument is almost independent of the core radius 
of each galaxy and the anisotropy parameter f3. We consider below the velocity dispersions 
at two different R values as a consistency check of our analysis. First however, we comment 
on another important feature of these figures which is relevant to a determination of galaxy 
parameters for lensing purposes. 

Kochanek (1996) has pointed out that the inclusion of a core radius will in general require, 
for self-consistency, a larger value of o"dm to be used in models, and has argued that this will 
cancel out the effect that non-zero core radii can have on suppressing lensing optical depths 
if the value of <tdm were instead kept fixed. Examining Figures 3 and 4 at approximately 0.5 
R e , where the different curves intersect, indeed illustrates the effect described by Kochanek. 
As the core radius is increased o"dm increases. However the key point is that the change 
in cdm is small compared to the huge variances in o"dm one sees in models with different 
values of /3, and also is small compared to the possibility of overestimating <tdm if central 
line-of-sight velocity dispersions are used. 

Note that in this section we chose 42 galaxies from 12 different data sets to be our galaxy 
sample, because for these we could find or estimate the velocity dispersions at both R c and 
0.54i? e • Among those 12 sets of observations, most use a slit size from one to three arcseconds. 
This should therefore not be significant when R e is on the order of ten arcseconds, as it is in 
this paper (see Table 1). Furthermore, from Figures 3 and 4, one can see that the velocity 
curves are relatively smooth between R c and 0.54i? e (i.e., theoretically, o"dm = 1-31o\r c = 
1.25cro.54R e ). This indicates that even if the slit size were up to half R c , this effect might still 
not strongly bias the results. 

3. 1 Velocity dispersions at R = R c 

Based on values of R c from RC3, we can estimate velocity dispersions <tdm of these 42 
galaxies. We also consider R e from various different papers and translate the different values 
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of velocity dispersions into the uncertainty of the velocity dispersion measured at R e , as 
given by RC3. We would quote either this uncertainty or the uncertainty from the observa- 
tions, whichever is larger. The other problem is that different observers might use different 
methods to measure velocity dispersions. We calculate the weighted mean of all these ve- 
locity dispersions for one galaxy as our final line-of-sight velocity dispersion of that galaxy 
(Lauer 1985,Bevington & Robinson 1992). We list the information of i? c ,B T (0), distances 
and velocity dispersions of all 42 galaxies in Table 1. 

Using Table 1 and a suitable value for = -19.66 + 51og 10 h ± 0.30 as discussed in 

Section 2, we can calculate the least-square-fit of the relation between log 10 (L/L*) and 
log 10 (a Ro /km s" 1 ) (Figure 5). We exclude galaxies NGC 4251, NGC 7796, NGC 4486B, and 
NGC 6411, because the former two galaxies do not have available distances from 7s, and the 
latter two galaxies are extreme outlyers (One galaxy is the farthest right one and the other 
galaxy is the lowest one in Figure 5.), and there is no available observational value of R c of 
NGC 4486B in RC3. We obtain 

log 10 (L/L*) = (-4.04 ± 0.49) + (1.89 ± 0.22) log 10 K c /km s" 1 ) (13) 

with x 2 = 48 for 38 galaxies in Figure 5. We use both uncertainties of luminosities and 
velocity dispersions to work out uncertainties in equation (13) and its \ 2 ■ If we se t L = L* 
in equation (13), then we find the velocity dispersion of luminous elliptical galaxies at the 
effective radius, a* Rc = 136 ± 15 km s" 1 . From the discussion in Section 3, we then multiply 
this number by 1.31, to get a*^ M 1.31<t|j o ps 178 km s _1 . 

In order to consider the confidence levels of the velocity dispersion we determine here, we set 
Ax 2 = 3.5 for the 68 per cent confidence level of velocity dispersion, and 7.8 for the 95 per 
cent confidence level (Lampton, Margon & Bowyer 1976), because we actually have three 
adjustable parameters (cr* Rc , L*, and the power law between or c and L/L*), even though we 
only vary a* R and fix the other two parameters. The final 95 per cent confidence level is then 
determined to be <t^ m = 178^8 k m s -1 . 

Several issues are relevant to this result. In the first place, note that the L vs a relation in 
the above equation differs from the standard Faber- Jackson relation (Faber & Jackson 1976). 
Note however that this FJ relation is appropriate for central velocity dispersions. In addition, 
it may be there case that the elliptical galaxies do not form a uniform population, but rather 
that bright galaxies and faint galaxies follow separate Faber- Jackson curves (this point was 
raised to us by J. Peebles, who has been investigating this issue (Fukugita & Peebles 1999)). 
Nakamura and Suto (1997) have also pointed out how the value of the power law relation can 
differ from 4. It thus may not be appropriate to enforce this relation on the bulk sample. To 
explore this latter possibility, and because it is galaxies with the largest velocity dispersions 
which will dominate in the analysis of lensing statistics, we considered a subset of our galaxy 
sample with log 10 ((XR o /km s _1 ) > 2.2 (i.e., 23 galaxies) and rederived the L-a relation. In 
this case, we find: 

log 10 (L/L*) = (-6.20 ± 1.82) + (2.83 ± 0.79) log 10 K e /km s" 1 ) (14) 
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Fig. 5. Least-square-fit of log 10 (L/L*) vs. log 10 (cTR e /km s _1 ). The center solid line is the best fit, 
as shown in equation (13). The upper dash line is the best fit plus 0.09, and the lower dash line is 
the best fit minus 0.09. These two lines are the 95 per cent confidence limits. 

with x 2 — 23.6 and a* Rc = 155 km s -1 . In this case, we find a somewhat higher value of 
°dm = 203 km s _1 as might be expected given the features of this subsample. 

A recent result from galaxy-galaxy weak lensing (Fischer et al. 2000) constrains the velocity 
dispersion a v of the average foreground galaxy to be 150-190 km s _1 and this result translates 
to a odm (using our notation) of the average foreground galaxy to be 210-266 km s -1 . It 
is interesting to compare this result with the estimate we would derive here. Based on the 
luminosity of the average foreground galaxy given in Fischer et al. (2000), i.e., 8.7x 10 9 L SMn in 
the g' band, our formula would predict a <jdm around 220 km s -1 , which is in good agreement 
with the range of <7dm determined by Fischer. 



3.2 Velocity dispersions at R = 0.54R C 

Another fit between luminosities and velocity dispersions of sample galaxies has been calcu- 
lated at R = 0.5AR C . Velocity dispersions at R — 0.54i? e are estimated and listed in Table 1. 
We still exclude NGC 4251, NGC 7796, and NGC 4486B, for the same reasons described in 
the above section. We find 

log 10 (L/L*) = (-5.01 ± 0.55) + (2.28 ± 0.24) log 10 K 5 4/? c /km s" 1 ) (15) 
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Fig. 6. Least-square-fit of log 10 (L/L*) vs. log 10 (<7o.54,R e /km s _1 ). The solid line is the best fit, as 
shown in equation (15). The upper dash line is the best fit plus 0.08, and the lower dash line is the 
best fit minus 0.08. These two lines are the 95 per cent confidence limits. 

with x 2 — 47 and o"o. 54 ^ c = 159 ± 15 km s" 1 (Figure 6). As we have discussed in Section 3, 
if we multiply this number by 1.25, then we should get the expected cx5 M . Therefore, a^ M ~ 
l-25(7Q 54iJe = 198 km s" 1 in this case. If we consider a subset of our galaxy sample with 
log 10 (o r o.54_R c /km s _1 ) > 2.2 (i.e., 28 galaxies), we find 

log 10 (L/L*) = (-6.43 ± 1.34) + (2.88 ± 0.58) log^aoWkm s" 1 ) (16) 

with x 2 — 23 and 0Q_ 54Re = 169 km s _1 . In this case, we find a^ M = 211 km s -1 . 

We see that the two slopes in L-a relation obtained at different values of R e are close to 
each other. Furthermore, the inferred value of a^ M at it! = 0.54i? e is somewhat larger than 
the value of <rJj M at R = R c . This can be explained as follows. The above estimates are 
for the case of a purely finite isothermal distribution. If one adds to this distribution some 
central mass, such as a large central black hole, this will further increase the central velocity 
dispersion, and also change the relationship between <7dm and the line-of-sight velocity dis- 
persion at R c and half R e . By measuring the velocity dispersions at both points, however, 
one can hope to extract out the central mass contribution and also the isothermal contribu- 
tion. (Alternately, it is clear that the velocity dispersion at R c will be less sensitive to the 
former contribution, and thus can be used to approximate the isothermal contribution.) It is 
important to recognize that the central mass contribution can affect the velocity dispersion 
more significantly than it affects the optical depth for lensing, so that incorporating this 
contribution in an effective isothermal sphere velocity dispersion will lead to an overestimate 
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of the optical depth. 



As a quick demonstration, let us see how we can extract out the central mass contribution. If 
we simply assume v{x) oc l/x 4 due to a point mass at the center of galaxy, from equations (11) 
and (12), we can obtain 



where M is the mass of that point mass. Then, we can set up equations as follows: 



Plosive) - V\ OS ,g( R e) + g^T 

0^(0.54^) = < g (0.54it! e ) + 6 - ™ ARc) (18) 

where (Ji OSig (i?) is the isothermal contribution to the velocity dispersion of that galaxy as a 
function of R. As discussed above, we know that 1.31<Ji OSig (i? c ) = 1.25(7i OS)g (0.54i? e ) = cxdm- 
Thus, the appropriate isothermal effective <tdm can be calculated for each galaxy and the 
value of ctJ, m can be found through a least-square-fit method. A very quick estimation of the 
appropiate c^m from equations (14) and (16), is ct^m ~ 190 km s _1 , with the slope of the 
L vs. cr DM relation being 2.8. It is interesting to note that when values close to these are 
utilized in an analysis of lensing statistics (Cheng & Krauss 2000) one obtains mean angular 
splittings that are comparable to the observed splittings. 



4 DISCUSSIONS 



First note that both M^^j and a, fitted to the Schechter function are close to the results in 
Loveday et al. (1992). This is expected, because the B T and Bj bands are close. In fact, we 
would expect our Mg T , j to be slightly fainter than (= -19.71 + 51og 10 /i), because the 
Bj band is measured at a longer wavelength and most galaxy luminosities are peaker toward 
the redder end of the spectrum. On the other hand, our Mg T ^ ^ is brighter than Mg T (= 
-19.37 + 51og 10 h) given in EEP. Given the Galactic extinction correction, this is also what 
we would expect. We have also shown that the Schechter a of elliptical galaxies is shallower 
than that for total galaxies, which is consistent with EEP and Loveday et al. (1992). 

The ratio of our 0* to the one in Loveday et al. (1992) is about 9.6 per cent. In comparison 
with the 0* in Mobasher, Sharpies & Ellis (1993), the ratio is about 12 per cent. These 
two numbers indicate the fraction of elliptical galaxies out of total galaxies, and they are 
consistent with the conclusion of Postman & Geller (1984), of 12 ± 2 per cent for this ratio. 
Our results here suggest that the population of elliptical galaxies is (10.8 ± 1.2) per cent of 
all galaxies. 
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Regarding our luminosity function determination, it is of course still possible that the galaxies 
in our sample set do not form a complete sample set. Currently, we calculate < > = 
0.45 (Schmidt 1968). If we had an incomplete elliptical galaxy sample set, and if we would 
believe those 'missing' E galaxies are near apparent magnitude 12.19, then we would expect 
< tj^— > larger than 0.45. Also, according to the theoretical prediction in Zucca, Pozzetti & 

V max 

Zamorani (1994), we would then have an even a dimmer M* and a steeper a. A fainter M* 
would result an even smaller a^, M than the value calculated in Section 3.1 and Section 3.2. 

The analysis in Section 3.1, suggests o* Rc = 136± 15 km s -1 . It implies that <rfj M = 178± 
29 km s _1 . This uncertainty is large, because the uncertainties from distances of galaxies in 
the sample set and from M^™ are large. (Also note that the best fit value if we account 
for possible central mass concentrations is u^ M ~ 190 km s -1 , which falls within the range 
included above.) We also notice that M^m is a sensitively coupled to o"Jj M . For example, 
if we use Mg T(0) = -19.9 + 51og 10 h (Fukugita & Turner 1991,Kochanek 1994), then we will 
increase a^ M by 8%-12%. On the other hand, if we use Mg T ( ) = -19.5 + 5 log 10 h, then we will 
decrease a^M by 5%-8%. One should be careful about this issue in any strong gravitational 
lensing analysis. The completeness of our galaxy sample discussed in Section 2.1 does not 
seem to lead to a significant uncertainty by comparison. For example, different numbers of 
galaxies included in our sample produce a variation of 4%-6% in (j DM , which is within the 
values estimated from the uncertainty of M* itself. It might be interesting to consider the 
relation between log 10 (L/L*) and log 10 (o"/km s^ 1 ) in other bands (e.g. V or R band, if the 
observation is possible; or a radio survey), to see whether including more elliptical galaxies 
would reduce the intrinsic scatter of this power law relation. If so, then the improved value of 
cx^M an d the power law from that survey would be a great help in the study of gravitational 
lensing statistics. Based on our present result, the optical depth of gravitational lenses would 
be expected to be reduced by 43 per cent compared to the one in Kochanek (1993), if all 
other parameters are unchanged. 

In this paper, we treat the anisotropy parameter (3 as a constant. We expect that our con- 
clusion will not change too much even if (5 is a function of radial distance r, as long as < 
0.4 (Kochanek 1994) (Figure 3 and Figure 4). 

Finally, the theoretical modelling in Section 3, suggests that of os /o"DM * s strongly dependent 
on models and parameters within 0.1i? e (Figure 3 and Figure 4). We also observe that the 
ratio (J\ os /cTi)M converges at about 0.5i? e for all different parameters. This suggests that we 
should use observational values of <J\ os /<Ji)M between 0.4i? e and R e , if we want to reduce 
the intrinsic scatter from theoretical modelling. We hope observers will consistently measure 
velocity dispersions of elliptical galaxies at both R c and 0.54i? e in order to refine our results. 
As we have discussed in Section 3.2, one can extract extra velocity dispersions at R c (and 
0.54-Re), and get an effective value of <tdm for each galaxy. As data improves, we expect 
that the techniques described here will be useful in further constraining input parameters in 
strong gravitational lensing analysis. 
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Name 


Rc(") 


Source 


m(B T (0)) 


Distance (/i -1 Mpc) 


<?R e (km s 1 ) 


o"0.54fl o (km s~ 


NGC 636 


19.4 


FIH 


12.16 


15.0 ± 1.8 


123 ± 10 


128 ± 3.4 


NGC 720 


36.1 


BDI 


11.13 


20.5 ± 4.4 


174 ± 12 


197 ± 4.7 


NGC 777 


34.4 


JS 


12.24 


44.0 ± 3.8 


275 ± 92 


228 ± 16 


NGC 1052 


33.7 


BDI/DB/FI 


11.33 


17.2 ± 3.7 


176 ± 3.4 


182 ± 2.1 


NGC 1399 


40.7 


SBS/WF 


10.22 


14.22 ± 0.88 


230 ±11 


230 ± 20 


NGC 1400 


29.3 


BZ 


11.87 


19.9 ± 1.9 


170 ± 25 


180 ± 6.0 


NGC 1404 


23.8 


FIH/WF 


10.88 


14.22 ± 0.88 


209 ±8.1 


229 ± 3.6 


NGC 1700 


18.5 


BSG/FIH 


12.00 


31.4 ± 4.7 


195 ± 7.9 


221 ± 6.4 


NGC 2434 


31.4 


CQ 


11.61 


19.7 ± 4.2 


190 ± 14 


190 ± 15 


NGC 2513 


32.9 


JS 


12.48 


46.6 ± 9.9 


185 ± 82 


221 ± 15 


NGC 2663 


54.6 


CQ 


10.23 


22.4 ± 4.7 


220 ± 14 


257 ± 17 


NGC 2778 


15.7 


JS 


13.21 


36.3 ± 7.7 


71 ± 26 


88.7 ± 8.1 


NGC 3115 


32.1 


BSG 


9.75 


10.2 ± 2.2 


100 ± 10 


110 ± 15 


NGC 3193 


26.7 


BSG 


11.73 


24.6 ± 3.7 


145 ± 25 


200 ± 20 


NGC 3377 


34.4 


BSG/FI 


10.98 


8.6 ± 1.3 


79.2 ± 4.8 


106 ± 7.5 


NGC 3379 


35.2 


DB/SBS 


10.17 


8.6 ± 1.3 


186 ± 6.7 


197 ± 4.8 


NGC 3557 


30.0 


FI 


10.79 


24.0 ± 5.1 


208 ± 15 


216 ± 8.5 


NGC 3608 


33.7 


JS 


11.71 


19.9 ± 2.4 


165 ± 26 


143 ± 7.3 


NGC 3610 


15.4 


BSG 


11.58 


21.3 ± 4.5 


129 ± 11 


147 ± 5.7 


NGC 3706 


20.3 


CQ 


11.87 


30.4 ± 6.4 


238 ± 8.9 


257 ± 8.3 


NGC 4251 


17.7 


BSG 


11.54 


(10.14 ± 0.56) 6 


85 ± 10 


110 ± 15 


NGC 4261 


36.1 


DB 


11.32 


27.8 ± 5.9 


248 ± 5.9 


292 ± 4.4 


NGC 4278 


34.4 


DB 


10.96 


14.7 ± 2.2 


178 ± 12 


197 ± 4.3 


NGC 4291 


17.3 


BSG/JS 


12.32 


30.3 ± 3.7 


187 ± 18 


200 ± 6.5 


NGC 4472 


104 


SBS/SZ/WF 


9.26 


13.33 ± 0.71 


198 ± 9.2 


261 ± 14 


NGC 4486B 


(3.07) a 


BSG 


14.26 


13.33 ± 0.71 


167 ± 8.3 


190 ±7.1 


NGC 4486 


94.9 


WF 


9.47 


13.33 ± 0.71 


238 ± 12 


(265 ± 2.6) d 


NGC 4494 


48.7 


JS 


10.61 


7.0 ± 1.5 


127 ± 21 


135 ± 8.2 


NGC 4564 


19.8 


BSG 


11.91 


13.33 ± 0.71 


80 ± 10 


105 ± 10 


NGC 4621 


40.5 


Bender 


10.5 


13.33 ± 0.71 


140 ± 40 


210 ± 10 


NGC 4660 


12.2 


BSG 


12.12 


13.33 ± 0.71 


110 ± 10 


135 ± 7.1 


NGC 5018 


22.8 


CQ 


11.34 


29.8 ± 3.1 


210 ± 7.1 


215 ± 14 


NGC 5576 


18.1 


BSG 


11.71 


16.5 ± 2.5 


100 ± 20 


185 ± 10 


NGC 5812 


25.5 


BZ 


11.87 


21.1 ± 4.5 


173 ± 117 


175 ± 17 


NGC 6411 


28.6 


JS 


12.64 


27.8 ± 5.9 


333 ± 120 


167 ± 15 


NGC 7144 


32.1 


SZ 


11.68 


18.4 ± 2.7 


155 ± 14 


145 ± 15 


NGC 7454 


25.0 


BSG 


12.7 


16.2 ± 3.4 


97.7 ± 8.3 


130 ± 10 


NGC 7507 


30.7 


BZ 


11.28 


17.5 ± 3.7 


180 ± 15 


207 ± 10 


NGC 7626 


38.6 


JS 


12.06 


35.8 ± 3.8 


166 ± 86 


208 ± 17 


NGC 7785 


23.3 


BSG 


12.41 


45.0 ± 9.6 


210 ± 30 


192 ± 8.9 


NGC 7796 


21.2 


BZ 


12.39 


(32.9 ±0.24) i > 


225 ± 37 


233 ± 12 


IC 179 


16.5 


JS 


12.98 


44.0 ± 3.8 


124 ± 21 


171 ± 8.7 



a This value is from BSG. 

6 The local peculiar velocity contributes to the distance estimate in RC3. 
c Reference WF (Winsall & Freeman 1993) is not used. 
d This value is from DB. 

Table 1 

The galaxy sample used in deriving the velocity dispersion-luminosity relation. The names of galax- 
ies are listed in the first column. Their corresponding effective radii, observational sources, apparent 
magnitudes, distances, velocity dispersions at effective radii, and velocity dispersions at 0.54 effec- 
tive radii are listed in columns 2 to 7, respectively. 



